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Abstract 

An algorithm for quantum computing Hamiltonian cycles of simple, 
cubic, bipartite graphs is discussed. It is shown that it is possible to 
evolve a quantum computer into an entanglement of states which map 
onto the set of all possible paths initiating from a chosen vertex, and 
furthermore to subsequently project out all states not corresponding to 
Hamiltonian cycles. 

A Hamiltonian cycle is a path on a graph which visits each vertx l..n exactly 
once, returning to the original vertex in the final step. Here we will discuss only 
simple (without loops or multiple edges), cubic (each vertex has only 3 edges), 
bipartite (each black vertex is connected only to white vertices) planar graphs. 
Such a graph is shown in Figure 1(a), with a Hamiltonian cycle indicated by 
the darker edges. 

Classically no efficient algorithm exists to resolve the question of whether 
any given such graph has a Hamiltonian cycle, although over the years many 
results in Graph Theory have isolated certain special cases. The fundamentals 
of Graph Theory and many associated classical algorithms are well explained 
in jjj, and a comprehensive exposition of Hamiltonian cycles and the related 
Travelling Salesman Problem can be found in Q . A good summary of quantum 
computing can be found in || . 

We shall see that the added power afforded us by a quantum computer's 
ability to carry through parallel computations in a single step enables us to 
compute all possible paths on a given graph. To achieve this we require n 
registers each composed of n qubits, each qubit corresponding to one vertex. 
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Three qubits of the first register, denoted a, are involved in every step of the 
algorithm. The other n — 1 registers contain qubits which in the main are in the 
state |0), except for the qubit corresponding to the walkers current position. 

To clarify the concepts consider the trivial but illustrative case of the square 
in Figure 1(b). We imagine a walker starting at vertex 1 and so prepare our 
quantum computer with register a in the state |1, 0, 0, 0) a , register 1 in the state 
1 1, 0, 0, 0) i = |l)i and the remaining registers empty. We envisage a controlled 
unitary operation R? for the jth step from vertex i. i?| is conditioned on the 
ith qubit in register j being in the state 1. It acts onthe 2 qubits in registers a 
and j + 1 which correspond to vertices adjacent to i. In a series of three steps 
we wish the initial state |1, 0, 0, 0) a |l)i to evolve as follows (empty registers not 
shown and register labelling dropped), 

R\ : {|1,1,0,0)|1}|2>; |1, 0, 1, 0)|1)|3)} 

R\R\: {|0,1,0,0)|1)|2)|1); |1, 1, 1, 0)|1)|2)|3); |0, 0, 0, 1)|1)|4) |1); 

|1,0,1,1)|1)|4)|3)} 

R\R\: {|0,0,0,0}|1}|2)|1)|2>; |0, 1,0, 1)|1)|2)|1)|4); |1, 0, 1, 0)|1)|2)|3)|2); 

|0,0,0,0)|1)|4)|1)|4); |0,1,0,1)|1)|4)|1)|2); |1, 0, 1, 0)|1)|4)|3)|4); 
|1,1,1,1)|1)|2)|3)|4); |1,1,1,1)|1)|4)|3)|2)}. 
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Figure 2: A schematic showing a suitable form of operator Uf , where a, b, c label 
the three vertices adjacent to i in the graph being considered. 



The {}'s denote a superposition of those states enclosed (ignoring normal- 
isation and phases for the moment). We see that only those paths which are 
possibly Hamiltonian cycles in the next step contain all l's in register a. This 
is because stepping back to an already passed vertex changes a pre-existing 1 
at that site to a 0. Thus projecting out this state of register a will leave the 
quantum computer in an entanglement of states which correspond to the graph's 
Hamiltonian cycles, if they exist. 

To progress to the slightly harder problem of cubic graphs we need to be 
more specific about the form of the unitary transform Uf required. It will act on 
6 qubits, and also be conditioned on qubit i in register j being 1. It will apply 
elementary NOT type operations to the 3 qubits a, 6, c in register a which label 
the vertices adjacent to i, and write l's into the same sites in the (previously 
empty) register j + 1 similarly to that discussed above. A schematic of such 
a transform is shown in Figure 2, using a version of Feynman's [Q notation 
developed in ||. The © is the elementary (1 bit) NOT operation given by the 
matrix Time progresses from left to right. The 3 qubit transformation V 

is given in the lexicographically ordered basis |0,0, 0);|0,0,1),...|1,1,1) by the 
matrix 
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As an example of how the computation would proceed, consider again Fig- 
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ure 1(a). A first application of U\ would evolve the initial state |1,0..0) into 
the superposition ^=(|1, 1, 0..0)|1) |2) + |1, 0, 0, 1, ..0) |1) |4) + |1, 0..1) |1) |8)). Sub- 
sequently we would apply R^R\R\ = R\ 4 8 followed by R\ 3 5 7 , then R\ 4 6 g 
and so on until we finally apply R 7 1, 3, 5, 7. Projecting out those states which 
contain all l's in register a will leave the computer in a superposition of states 
which map onto Hamiltonian cycles. A final measurement will reveal one of 
those cycles, should they exist for the particular graph under consideration. 

It is easy to see that in general we require 0(n 2 ) applications of Uf . This 
is obviously impractical with present technology. However algorithms such as 
the one above help us understand what quantum computers are capable of in 
principle and are therefore important in understanding the relation between 
quantum and classical complexity classes. 
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